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We present a numerical study of the magnetization process of frustrated quantum spin-5 chains 
with 5=1, 3/2, 2 as well as the classical limit. Using the exact diagonalization and density-matrix 
renormalization techniques, we provide evidence that a plateau at one third of the saturation mag- 
netization exists in the magnetization curve of frustrated spin-5 chains with S > 1/2. Similar to the 
case of 5=1/2, this plateau state breaks the translational symmetry of the Hamiltonian and realizes 
an up-up-down pattern in the spin component parallel to the external field. Our study further shows 
that this plateau exists both in the cases of an isotropic exchange and in the easy-axis regime for 
spin-5=l, 3/2, and 2, but is absent in classical frustrated spin chains with isotropic interactions. 
We discuss the magnetic phase diagram of frustrated spin-1 and spin-3/2 chains as well as other 
emergent features of the magnetization process such as kink singularities, jumps, and even-odd ef- 
fects. A quantitative comparison of the one-third plateau in the easy-axis regime between spin-1 
and spin-3/2 chains on the one hand and the classical frustrated chain on the other hand indicates 
that the critical frustration and the phase boundaries of this state rapidly approach the classical 
result as the spin 5 increases. 



I. INTRODUCTION 

The search for novel quantum phases has been a strong 
motivation for theoretical investigations of quantum spin 
models. Quantum states that do not resemble classically 
ordered magnetic states may be realized as a consequence 
of reduced dimensionality or competing interactions. Fa- 
mous examples are the Haldane phase of integer spin 
chains, 1 or the spin liquid ground states of spin- 1/2 lad- 
der systems, 2 both of which exhibit a finite spin gap. 3-5 
The presence of competing interactions causes frustration 
and drives quantum phase transition, 6 as for instance in 
the frustrated spin-1/2 chain that undergoes a transition 
from a critical, gapless phase into a spontaneously dimer- 
ized one. 7,8 Of great interest are the excitations above 
these ground states, which are reflected in experimentally 
accessible quantities such as thermodynamic properties. 

Here we are concerned with the interplay of frustra- 
tion and an additional external parameter, the magnetic 
field. A magnetic field can close zero-field gaps, but 
may as well induce gapped magnetic excitations in fi- 
nite fields. 9,10 The magnetic phase diagram can exper- 
imentally be mapped out by studying the magnetiza- 
tion process, which theoretically may be computed by 
means of powerful and flexible numerical methods such 
as the density matrix renormalization group technique 
(DMRG). 11,12 One of the intriguing features of the mag- 
netization curve of low-dimensional quantum magnets is 
the emergence of plateaux at finite magnetic field, which 
indicate the presence of massive excitations. This phe- 
nomenon has been predicted for several low-dimensional 
spin models. 9,10,13-29 Experimentally, it has been ob- 
served in many quasi-low dimensional magnetic mate- 
rials, realizing networks of spin-1/2 (Rcfs. 30-34) as well 
as of spin-1 moments. 35 It is the purpose of this work to 



study the magnetization process of frustrated quantum 
spin chains with 5 > 1/2 and establish the existence of 
a plateau state at one third of the full magnetization M. 

The possible existence of magnetization plateaux in 
quantum spin chains was predicted by Oshikawa, Ya- 
manaka, and Affleck. 9 A notable result of this work pro- 
vides a necessary condition for the existence of plateau 
states at finite M: 

P S(i-M)eZ. (l) 

Here, 5 is the spin, p denotes the period (or length of 
the elementary unit cell) of the plateau state in real 
space. The theorem implies that M always has a rational 
value on any plateau and it allows for translational in- 
variance to be broken spontaneously, i.e., p can be larger 
than 1. Note that unfrustrated spin-5 chains usually 
show a smooth magnetization curve, as has been known 
for quite some time already, 36 implying that compet- 
ing interactions beyond a simple nearest neighbor model, 
anisotropics, or special geometries are responsible for the 
plateau formation. Indeed, theoretically, the existence 
of plateaux has been established for, e.g., 5=3/2 chains 
with onsite anisotropy 9,16,21 at M=l/2, in frustrated and 
dimerized spin-1/2 chains at M=l/4, 13 three leg ladders 
at M=l/3, 10 or exotic models such as the orthogonal 
dimer chain. 23 

The simplest model with competing interactions in one 
dimension is a spin-1/2 chain with nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) exchange interactions, 
also called a zigzag ladder. Only quite recently the ex- 
istence of a M=l/3 plateau in the magnetization curve 
of frustrated antiferromagnetic 5=1/2 chains has been 
revealed by DMRG calculations. 2 The plateau is accom- 
panied by broken translational symmetry with a period 
p = 3 and an up-up-down (uud) structure in the spin 
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component parallel to the external field. While several 
higher-5 chain models have been studied in the context 
of magnetization plateaux, 9. 3-5, 16,20,21,37 ^he case of frus- 
trated chains with 5 > 1/2 has not been addressed in suf- 
ficient detail. In particular, the existence of the M=l/3 
plateau in these systems, which we report in this work, 
has not been explored. Our main result, based on exact 
diagonalization (ED), DMRG, classical Monte-Carlo sim- 
ulations (MC), and spin- wave theory, is that the M=l/3 
plateau is realized in frustrated quantum spin chains with 
5 > 1/2 in the case of isotropic exchange and in the easy- 
axis regime as we show for the specific examples of 5=1, 
3/2, and 2. As for spin-1/2, 24 this plateau always has 
broken translational symmetry with period p = 3 and 
the uud structure. 

A M=l/3 plateau with this Neel-type of order emerges 
in the Ising- limit of both quantum 37-39 and classical frus- 
trated chains (see Ref. 40 anf further references therein). 
While the classical frustrated chain does not exhibit a 
finite magnetization plateau in the case of isotropic in- 
teractions as we show in this work, the M=l/3 plateau 
is stable against quantum fluctuations. The parameter 
region, however, where such plateau exists for the case of 
isotropic interactions shrinks as the spin 5 grows. As we 
compute the full magnetization curves, we also obtain in- 
formation on other intriguing features such as kink singu- 
larities, macroscopic jumps, and even-odd effects in the 
magnetization curve, and a very rich phase diagram is 
indeed found. We present a qualitative discussion of the 
magnetic phase diagram of spin-1 and spin-3/2 chains, 
which extends the picture emerging from previous stud- 
ies for zero 41-46 and finite magnetic fields. 47 

In contrast to the case of isotropic interactions, we find 
that the M=l/3 plateau region prevails in a substan- 
tial part of the magnetic phase diagram in the easy-axis 
regime. Its boundaries rapidly approach the classical re- 
sult in the easy-axis regime. To this end, we study the 
magnetic phase diagram of the classical frustrated chain 
as well by means of MC simulations and linear spin- wave 
theory, aiming at a comparison of emergent phases with 
the quantum cases of 5=1 and 5=3/2. 

In a very recent experimental study, a M=l/3 plateau 
has been observed in the frustrated diamond spin-1/2 
material Cu3(C03)2(OH) 2 (Ref. 33). The existence of a 
1/3 plateau has also been reported for the spin-1/2 trimer 
compound Cu3(P20eOH) 2 (Ref. 34). Our results may 
be of particular relevance for materials that have been 
suggested to realize frustrated spin-1 chains. One ex- 
ample is CaV205 with competing next and next-nearest 
neighbor antifcrromagnctic interactions. 48 A promising 
family of materials realizing higher spin-5 zigzag chains 
is NaR(WO<4)2 where R represents In, Sc (5=0), Cr 
(5=1/2), and V (5=0 or 5=1). 49 The exchange con- 
stants in the latter case are estimated to be of the or- 
der of 180K, 49 which may be small enough to access a 
substantial part of the magnetization curve. Finally, the 
formation of spin-1 zigzag chains has recently been re- 
ported for TI2RU2O7, for which a Haldanc gap of 110K 



is found. 50 We hope that our results will stimulate the 
search for the one-third plateau in potential frustrated 
spin chain materials. 

The rest of the paper is organized as follows. In Sec. II, 
we introduce the model and briefly discuss the numerical 
techniques employed in this work: DMRG, Lanczos, and 
MC simulations. In Sec. Ill, we present our numerical 
results for the magnetization process of frustrated spin 
chains with 5=1, 3/2, and 2 and isotropic exchange in- 
teractions. We demonstrate the existence of the M=l/3 
plateau for these quantum spins and numerically deter- 
mine the phase boundaries of the plateau region for spin 
1 and 3/2 in the field vs frustration plane. We study 
the magnetization process of spin-1 chains in more de- 
tail, focusing on (i) kink singularities in the magnetiza- 
tion curve, (ii) low-field phases, and (iii) and the com- 
bined effect of frustration and onsite anisotropics. The 
latter is important as most materials that realize Hal- 
dane chains have a substantial single-ion anisotropy. In 
Sec. IV, we turn to the case of anisotropic frustrated 
chains in the easy-axis region. The main purpose of this 
section is to demonstrate how the M=l/3 plateau ob- 
served in frustrated quantum spin-5 chains is connected 
to the classical limit. To this end, we first analyze a clas- 
sical frustrated spin chain by means of MC simulations, 
which allows us to map out the ground-state phase di- 
agram of the classical model in a finite magnetic field. 
Analytical results for several phase boundaries are pre- 
sented. We then perform DMRG calculations for frus- 
trated spin-1 and spin-3/2 chains at (i) a fixed exchange 
anisotropy, varying the frustration, and (ii) at a fixed 
frustration, varying the exchange anisotropy, and deter- 
mine the phase boundaries of the M=l/3 state for these 
two cases. A summary and conclusions are presented in 
Sec. IV. An Appendix is devoted to the details of a linear 
spin- wave calculation for the M=l/3 plateau state in the 
classical limit. 



II. MODEL AND METHODS 

The Hamiltonian of a frustrated spin-5 chain with on- 
site anisotropy D in a magnetic field h reads: 

H = Y,[Ji ^ ■ S i+ i + J 2 St ■ S J+2 - hSf + D(5?) 2 ] , (2) 

i 

where J\ > and J2 > are the antiferromagnetic NN 
and NNN exchange integrals, respectively Here, S, = 
(Sf,Sf,Sf) denotes a spin-5 operator acting on site i. 
We also include an exchange anisotropy A, the same for 
the NN and NNN interactions. Hence, in our notation, 
with the usual definition of raising and lowering operators 

S< ■ Sj = ±(S+Sr + Sr S +) + A 5f 5| . (3) 

We restrict our analysis to A > 1, and we set J\ = 1 and 
D = unless stated otherwise. The following normal- 



3 



ization for the magnetization is used [this also applies to 
Eq. (1)]: 



M=^f; St otal =^Sf, (4) 



where N is the total number of sites in the chain. Thus, 
717=1 at the saturation field h sa t. 

We compute the ground-state energies Eo(S* otal ,h = 
0) in subspaces labeled by S^ otal on chains with periodic 
boundary conditions (PBC) using the Lanczos algorithm. 
The ground-state energies of substantially larger chains 
with open boundary conditions (OBC) are calculated 
with DMRG. Typically we keep up to m — 400 states 
in our DMRG calculations. Then, we include the Zee- 
man term and obtain the field-dependent ground-state 
energies 

E a (S? otal , h) = E (Sl tal , h = 0)-h S* otal . (5) 

The magnetization curves are constructed by solving 
the equations E (S* otal ,h stcp ) = E (S% otal + s,h step ), 
which define those magnetic fields at which the mag- 
netization increases from M = S^ otal /(N S) to M' = 
(S* otai + s)/(N S). Steps larger than ,s=l may occur. 

In the classical limit, we treat Si as three-dimensional 
unit vectors given by the polar angle 0^ and azimuthal 
angle fa. For each value of h, we find the ground state 
of H using the "simulated annealing" classical MC al- 
gorithm. Starting from a random configuration of an- 
gles {Oi, fa} and an initial inverse temperature (3 = /3j n , 
we proceed as follows. We perform sweeps through the 
chain, and for each site i = 1, . . . N, two steps are per- 
formed: (i) The energy E = E({<di, fa}) of the current 
configuration is calculated from Eq. (2). (ii) New angles 
O'i = Oi+r-w and fa = (fi + r-w are tried, where r G [0, 1] 
is a random number and w is the maximal change al- 
lowed. The new angles are accepted with probability 
p = min{exp[/3(E - E')], 1}, where E' = is 
the energy of the configuration with 0^ and fa replaced 
by 0- and </>-. After each sweep through the chain, (3 is 
multiplied by the factor (^fl/An) 1 / , where /3fi is the final 
inverse temperature and G is the number of sweeps, w is 
adjusted such that the acceptance ratio is kept close to 
50% for efficiency. Typically, we use (3- m = 0.5, (3s = 10, 
and G = 10 5 . This slow-cooling technique prevents the 
simulation from getting stuck in local minima in the en- 
ergy landscape in most cases. Sometimes, close to phase 
boundaries, the system gets trapped in local minima, but 
repeating the algorithm a few times for the same set of 
parameters usually solves the problem. After /3g has been 
reached, we fine-tune the ground-state energy by repeat- 
ing steps (i) and (ii) with f3 — oo for several sweeps until 
the energy change after a sweep becomes negligible. 
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FIG. 1: (Color online) Magnetization curve M(h) for frus- 
trated chains with an isotropic exchange (A = 1) for (a) 5*=1 
(J 2 = 0.4), (b) 5=3/2 (J 2 = 0.4), and (c) S=2 (J 2 = 0.36). 
Solid(dashed) lines are DMRG(ED) results. The horizontal 
dotted lines mark M=l/3 and M=l. 



III. MAGNETIZATION PROCESS OF 
FRUSTRATED QUANTUM SPIN CHAINS WITH 
ISOTROPIC EXCHANGE INTERACTIONS 

In this section we present results for the magnetization 
curves of frustrated spin-S 1 chains with A = 1. As a main 
result of this paper, we provide numerical evidence that 
the M=l/3 plateau, previously observed in the magncti- 
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FIG. 2: (Color online) Onsite spin density {Sf)/S on 
the M—l/3 plateau state of frustrated spin-5 chains with 
5=1,3/2, and 2 and isotropic exchange interactions (A = 1) 
at J 2 = 0.4 (5=1 and 3/2) and J 2 = 0.36 (5=2). We show 
DMRG results for N = 120 sites. 



zation curve of frustrated spin-1/2 chains, 24,25 exists for 
higher spins as well. The boundaries of the plateau phase 
strongly depend on 5, and the region in the /i-J 2 plane 
where the M=l/3 state is realized shrinks as a function 
of S. Our MC simulations in the classical limit show 
no evidence for a magnetization plateau in the isotropic 
case. 



Magnetization curves of frustrated 5=1, 3/2, 
and 2 chains 



The magnetization curves at an intermediate frustra- 
tion strength are presented in Fig. 1 for the cases of spin 
1 [Fig. 1(a), J 2 =0.4], spin 3/2 [Fig. 1 (b), J 2 =0.4], and 
spin 2 [Fig. 1(c), J2=0.36]. The magnetization plateau at 
M=l/3 is clearly observed in all three plots, and it can 
thus be considered a universal feature of frustrated quan- 
tum spin chains. In accordance with the criterion Eq. (1) 
for the plateau formation, the translational invariance of 
the Hamiltonian is broken on the plateau, leading to a 
state with period p = 3 and an uud pattern in the onsite 
spin density (Sf) parallel to the external field. The lo- 
cal spin densities corresponding to the plateau state are 
shown in Fig. 2. 52 As we use OBC suited for DMRG cal- 
culations, one of the three possible patterns (up-up-down, 
up-down-up, down-up-up) is selected as a unique ground 
state due to energetically preferred orientations of edge 
spins, pointing up in our case. We further observe devi- 
ations from the perfect up-down-up pattern, selected in 
our case, at the boundaries: edge-spin induced spin-spin 
correlations decay into the bulk. 

Let us briefly comment on additional features of the 
magnetization curves presented in Fig. 1. In the case 
of spin 1 shown in Fig. 1(a), we identify at least six 
regions, labeled (i)-(vi) in the viewgraph. First, the 
Haldane gap 1,4,42 shows up as a zero-field plateau. A 



1/2 
1 

3/2 
2 

oc 



0.56 
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0.72 



TABLE I: Isotropic exchange, A=l: Lower and upper critical 
frustrations J2,crit,i/u f° r the formation of the M=l/3 plateau 
for spin 1/2 (from Ref. 24), spin 1, 3/2, and 2. 



specialty of open spin-1 chains is the presence of ef- 
fective spin-1/2 edge spins which have been studied in 
detail both for unfrustrated 4 and frustrated chains. 42 ' 51 
These edge spins couple to the so-called Kennedy triplet, 
which in the thermodynamic limit is degenerate with 
the ground state. 53 The existence of the effective free 
spin-1/2 edge spins emerges naturally within the AKLT 
(Afflcck-Kennedy-Lieb-Tasaki) description of the ground 
state of a spin-1 chain. 4 ' 54 Within our numerical pre- 
cision, this triplet causes a zero-field magnetization of 
M = 1/N seen in the Haldane phase. Similar edge spin 
effects are present in all open spin-S' > 1/2 chains (see, 
e.g., Ref. 51 and references therein): spin-3/2 chains 
have spin-1/2 edge spins whereas spin-2 chains have 
spin-1 edge spins, which consistently result in spurious 
M = 1/N and M = 2/N zero-field magnetizations seen 
in Figs. 1 (b) and 1(c). 

At h ps 0.82, the frustrated spin-1 chain enters a sec- 
ond region where the magnetization initially increases 
with a steep slope but then follows M(h) = \/h — b with 
b ps 0.74 for h > 0.9. The M=l/3 plateau is realized 
for 1.74 < h < 1.99. Note, however, that the phase 
boundaries of the plateau state cannot directly be ob- 
tained from TV = 120 sites for this frustration parameter, 
as finite-size effects are still significant in this regime. 
Both below and above the plateau (indicated as region 
iv), the magnetization increases in steps corresponding to 
A5 t 2 otal = 2. This may indicate binding effects between 
the elementary excitations in finite magnetic fields. In 
the high-field region below the saturation field, we high- 
light the presence of a kink in the magnetization curve at 
h ps 3.89, which correspond to jumps in the differential 
susceptibility dM(h) / dh. This kink splits the high-field 
region into parts (v) and (vi). A second, less obvious kink 
may exist at smaller fields h ps 2.8. The kink singulari- 
ties are related to middle-field cusp singularities recently 
studied in the literature. 47 ' 55,56 We will discuss the kink 
singularities in more detail in Sec. Ill C. 

Apart from the M=l/3 plateau, the prominent fea- 
tures of the magnetization curve of a frustrated spin- 
3/2 chain with J 2 = 0.4 are the zero-field gap and sev- 
eral kinks. At zero magnetic field, the frustrated S=3/2 
chain is gapless for J 2 < 0.29 and in a critical phase 
with antiferromagnetic correlations, while a gap opens 
beyond this value. 51 This transition at zero field is of the 
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FIG. 3: (Color online) Magnetic phase diagram of frustrated 
5=1 and 3/2 chains with an isotropic exchange (A = 1) as 
derived from ED and DMRG calculations of magnetization 
curves. Squares denote boundaries of 5=1 phases, while stars 
are for S=3/2. An arrow at J2,t ~ 0.75 at low fields marks the 
first-order transition from the (gapped) Haldane phase found 
for J2 < 0.75 to the likewise gapped double Haldane phase 
for the frustrated spin-1 chain. 59 



Kosterlitz-Thouless type and is similar to the behavior 
of the frustrated 5=1/2 chain. 7 ' 8 At larger J 2 > 0.29, 
the system enters a dimerized phase, with a finite spin 
gap seen in Fig. 1 (b). For the case of h > 0, several 
kink singularities exist in the magnetization process of 
spin-3/2 chains both below and above the plateau, as is 
indicated in Fig. 1(b). The derivative of the ground-state 
energies Eo(S^ oted , h = 0) with respect to M exhibits dis- 
continuous changes of the slope where the kinks emerge 
in the magnetization curve, indicating at second-order 
transitions. Similar to the case of spin 1, the magnetiza- 
tion increases in steps corresponding to 5 t z otal = 3 = 25 
around the plateau. This hints at a very rich phase dia- 
gram in a magnetic field that deserves to be explored in 
more detail. 

The magnetization curve of frustrated spin-2 chains at 
J-2 = 0.36 exhibits a zero-field Haldane gap, 45 similar to 
the case of spin 1, kinks at high fields, and the M=l/3 
plateau. 



B. Magnetic phase diagram of frustrated spin-1 
and 3/2 chains with isotropic exchange interactions: 
The M=l/3 plateau region 

We next discuss how the boundaries of the plateau re- 
gion depend on spin 5. First, we present results for the 
critical frustration for the plateau formation in Table I. 
Second, the phase boundaries are determined as a func- 
tion of frustration J2 for spin 1 and spin 3/2 and these 
results are summarized in Fig. 3. Finally, we present a 
discussion of additional phases present in the case of spin 
5=1 in a separate section, Sec. IIIC. 



Since the plateau is induced by frustration, we expect 
it to vanish as J 2 — > 0, as well as for J 2 3> J\- In the 
latter case, the system separates into two decoupled un- 
frustrated NN chains. For the extreme quantum case of 
5=1/2, DMRG studies 24 ' 25 report the existence of the 
plateau for 0.56 < J2/J1 < 1.25. On the other hand, 
as mentioned earlier, our MC simulations for a classical 
frustrated spin chain indicate the absence of the M=l/3 
plateau in the isotropic case A = 1, thus we expect the re- 
gion that allows for the uud state to become the ground 
state at M=l/3 to shrink with increasing spin 5. We 
therefore determine the critical lower and upper frustra- 
tion J2. C rit,i/u(5) for the formation of the plateau state. 
To this end, we compute the onsite spin density (5f ) for 
several J2, searching for the emergence of the uud pat- 
tern. We also perform a finite-size scaling analysis of the 
plateau boundaries and compare chains with open and 
periodic boundary conditions. This procedure results in 
estimates of the critical frustrations for spin-1/2 consis- 
tent with Ref. 24. Note that the precise determination of 
the critical frustrations is somtimes difficult as the gap, 
i.e., the width of the plateau, is typically exponentially 
small close to the critical frustrations. Our results are 
intended to study trends and we estimate the critical 
frustrations listed in Table I and II to be correct with an 
error of ±0.02. 

The results for 5=1, 3/2, and 2 arc summarized in Ta- 
ble I, together with data from Ref. 24 for spin 1/2. The 
data collected in Table I renders 5=1/2 a special case 
with critical frustrations very different from 5 > 1/2. For 
5 > 1/2, the critical lower J2, C rit,i(5) remains approx- 
imately constant, while J2, C rit,u(5) — J2,crit,i(5) shows 
a tendency to decrease with increasing 5, melting the 
plateau region from the large J2 side. 

We next elaborate more on the phase boundaries h = 
h{ J2) of the plateau state in the h- J2 plane for 5=1 and 
3/2. Our results from DMRG calculations for several 
J2 are summarized in Fig. 3, which shows the magnetic 
phase diagram of frustrated 5=1 and 3/2 chains at A=l. 
Additional magnetization curves for 5=1 and J2=0.0, 
0.3, 0.7, and 0.8 are shown in Fig. 4 and in Fig. 5 for 
5=3/2 and J 2 =0.3, 0.5, 0.6, and 0.8. 

Close to the upper critical fields J2, C rit, u (5), the deter- 
mination of the boundaries is hampered by the existence 
of additional steps on the plateau. An example is shown 
in the inset of Fig. 4(c) for J2 = 0.7, the steps are indi- 
cated by the arrow in the plot. Such steps have also been 
seen in the case of frustrated 5=1/2 chains, 24 and can 
be traced back to the open boundary conditions. First, 
as we have verified for several cases, in which such steps 
are observed, we find that the step height scales down 
to zero with increasing system size. Consistent with the 
interpretation that the steps are boundary effects, they 
are not seen if periodic boundary conditions are imposed. 
Finally, we compare the onsite spin density (5f ) on the 
plateau to that of the step states, i.e., the magnetiza- 
tions just below and the two magnetizations just above 
the plateau. The spin densities are shown in Fig. 6 for 
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FIG. 4: (Color online) Magnetization curves for frustrated 
S=l chains with an isotropic exchange (A = 1) for (a) J 2 = 0, 
(b) J 2 = 0.3, (c) J 2 = 0.7, and (d) J 2 = 0.8. Full lines are 
DMRG results for N = 120 sites, dashed lines are ED results 
for N = 15. The inset in panel (c) is an enlarged view of the 
M=l/3 plateau. The squares indicate the boundaries of the 
plateau as seen on N = 15 sites with PBC. 



J 2 = 0.7 and N = 120 sites. On the plateau (5f otal = 40), 
we find the up-down-up pattern. For S* otsd = 39,41,42, 
i.e., just above and below the plateau, one-domain wall 
excitation sits in the middle of the chain. 25 

In the magnetic phase diagram Fig. 3, we scale the 
field with the value of 5. We see that while the up- 
per and lower critical frustration J-x,{i,u) are similar for 
5=1 and 3/2, the width and the critical field, i.e., the 
boundaries of the plateau region exhibit a significant de- 
pendence on spin 5. In particular, we emphasize the 
non-monotonic dependence of the plateau's width on 5 - 
it is quite narrow for 5=1/2, 24 broader for 5=1 and 3/2, 
and it disappears in the classical case 5 = 00. 



C. Magnetic phase diagram of the frustrated 
spin-1 chain: Additional phases 

Finally, we present a discussion of additional phases 
present in the case of spin 5=1. To our knowledge only 
the magnetization curve of the NN spin-1 chain has been 
computed (see, e.g., Ref. 57). Based on a field-theoretical 
analysis, a first suggestion for the magnetic phase dia- 
gram of frustrated 5=1 chains has been put forward in 
Ref. 47 for the low- and high-field parts. Our results 
basically confirm their picture and provide additional in- 
formation on the so far unexplored middle-field region 
where the M=l/3 plateau is found. 

Several studies of frustrated spin chains have focused 
on the zero-field phases, 42,44-47,51,58,59 with special in- 
terest in the emergence of chirality. As is well known, 
unfrustratcd 5=1 chains exhibit a Haldane gap at zero 
magnetic field, in contrast to half-integer spin chains 
that have gapless excitations. 1,5,11 As earlier studies 
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FIG. 5: (Color online) Magnetization curves for frustrated 
5=3/2 chains with an isotropic exchange (A = 1) for (a) 
J 2 = 0.3, (b) J 2 = 0.5, (c) J 2 = 0.6, and (d) J 2 = 0.8. Full 
lines are DMRG results for N = 60 sites, dashed lines are ED 
results for N = 12 sites with PBC. 



have shown, frustrated 5=1 chains have a gap for any 
j 2 42,60,61 r £) ae Haldane phase persists up to J 2 ~ 0.75, 
and is then followed by another gapped phase, the 
doublc-Haldane region. 42 This transition is a consequence 
of the frustration driving the system from the unfrus- 
trated one-chain (J 2 — > 0) to the two-chain limit ( J2=l). 

The existence of a spin gap gives rise to a zero-field 
plateau in the magnetization curve, which is evident from 
our data shown for 5=1 and J 2 =0.0, 0.3, 0.4, 0.7, 0.8 
[Figs. 4(a), 4(b), 1(a), 4(c), 4(d)]. The zero-field mag- 
netization of M = 5 t 2 otal = 1/N due to the Kennedy 
triplet is seen in the Haldane phase, but disappears at 
the transition to the Double-Haldane phase. 42 

We therefore estimate the spin gap from the field 
where the total 5* otal increases from 5* otal = 1(0) to 
S'totai = 2(1) in the Haldanc(Doublc-Haldane) phase. 
Our results for the spin gap arc included in the magnetic 
phase diagram in Fig. 3 (squares) and are in agreement 
with the data of Ref. 59. 

While the magnetization curve of both the unfrus- 
trated chain depicted in Fig. 4(a) and of chains with large 
J2 (not shown in the figures) are mostly featureless, a 
more interesting behavior arises in the range of parame- 
ters close to the plateau region. Around the plateau itself 
we find evidence for a narrow region where the magne- 



total 



tization increases in steps corresponding to A5 ( 
states with an odd 5^ otal never become ground states. 
Such a region exists in frustrated 5=1/2 chains as well, 
both in its anti- and ferromagnetic versions. 24,63 The 
corresponding state has been called the even-odd (EO) 
phase. 24 While the EO phase for 5=1/2 is predominantly 
realized at large J2, we observe the EO effects only in a 
narrow region around J 2 « 0.4 [see Fig. 1(a)]. Starting 
from the limit of two decoupled chains, the existence 
of the EO regions can be motivated in a simple pic- 
ture. When Ji = 0, spins always flip in pairs, one on 
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FIG. 6: (Color online) Onsite spin density (5f ) state for 5=1, 
,h = 0.7, A = 1, and N = 120 sites (DMRG). Panels (a)-(d): 
5f otal =39-42. 5 t z otal = 40 corresponds to M=l/3. 



each chain. The nontrivial question is whether a finite 
J\ preserves this effect leading to pair-binding of excita- 
tions. This is indeed the case for the frustrated 5=1/2 
chain. 24 ' 63 For spin 1, however, Kolezhuk and Vekua 47 
conclude that an EO phase is not possible in frustrated 
chains. Our DMRG results do not contradict this predic- 
tion as we find EO effects in a region that has not been 
studied in Ref. 47. 

To gain a qualitative understanding of the behavior 
in high magnetic fields, it is instructive to consider the 
dispersion of a single magnon, i.e, one spin flipped with 
respect to the fully polarized state. The dispersion e k , 
generalizing results given in Rcfs. 55 and 62, reads 

e k = 2S[J 1 cosk + J 2 cos2fc - A(J 1 + J 2 )] -h. (6) 

Here, and throughout, fc denotes the momentum. We also 
have included the exchange anisotropy A [see Eq. (3)] 
for the sake of generality. Equation 6 first allows us to 
derive the saturation field, which depends on S only by 
an overall prefactor: 62 



h sa ,t 



25[Ji - J 2 + A(Ji + J 2 )] for J 2 < -j, 



25 



8J 2 



J 2 + A(Ji + J 2 ) 



for J 2 >^.(7) 



The saturation field is plotted as a dashed line in Fig. 3. 
For A = 1 the saturation field is independent of J 2 for 
J 2 < 1/4. 

For J 2 > 0.25, a kink singularity splits off the sat- 
uration field and divides the high-field region into two 
phases. A second kink emerges at lower fields at J 2 « 0.4 
and then approaches the first one. We observe the first 
high-field kink for 0.25 < J 2 < 0.6, where it then merges 
with the second kink and eventually disappears, at least 
on the system sizes studied here. Just below the satu- 
ration field, where the one-magnon state is expected to 
give a correct description, we can understand its emer- 
gence along the lines of Ref. 55, where a similar effect has 
been observed for spin-1/2 chains. It turns out that for 
J 2 > 0.25 the global minimum of e k at fc = 7r splits into 
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FIG. 7: (Color online) Magnetization curves for the frustrated 
5=1 chain with isotropic exchange, A = 1 and Ji = 0.6, and 
a finite onsite anisotropy. (a) D = —0.5, (b) D = 0.1, (c) 
D = 0.2, (d) D = 0.5. DMRG results are for N = 60 sites 
(straight lines); ED (dashed lines) for N = 15 sites and PBC. 



two local minima at the position of the high-field kink 
singularity. We can further conclude that the high-field 
region between this kink and the saturation field is in- 
commensurate with cos fc = — Ji/4J 2 . This agrees with 
Ref. 47 where this part of the phase diagram has been 
suggested to exhibit chiral order. 

We emphasize that the interpretation of the kink singu- 
larity based on the one-magnon dispersion Eq. (6) is justi- 
fied only just below the saturation field. Indeed, even for 
the case of spin 1/2, the nature of the two phases left and 
right of the kink singularity is currently under controver- 
sial discussion. According to Refs. 24 and 55, it signals 
a transition from a one-component Tomonaga-Luttingcr 
liquid (TL1) to a two-component Tomonaga-Luttingcr 
liquid (TL2) state. This interpretation has recently been 
challenged by Kolezhuk and Vekua. 47 Using bosonization 
techniques, they conclude that the kink can be explained 
in terms of a renormalization of the magnetization due 
to an irrelevant operator in the incommensurate region 
neighboring the TL2 (or chiral phase) on the large J 2 
side. 

Further studies are necessary to relate the three fea- 
tures of the M(h) curves - the plateau, kink singularities, 
and EO effects - with each other and should involve a de- 
tailed study of excitations of 5=1 and 3/2 chains as well. 
Also, a numerical investigation of the high-field phases 
including the possible existence of chiral order is a neces- 
sary step towards a complete understanding of the phases 
5=1 chains in finite magnetic fields. 



D. Magnetization process of frustrated 5=1 chains 
with onsite anisotropy 

Most materials that realize 5=1 chains typically also 
have a significant onsite anisotropy D ^ of the form 
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FIG. 8: (Color online) Magnetization curve of a frustrated 
classical chain (A = 2) for J2 = 0.8 obtained by means of 
MC simulations (circles). The horizontal, dotted line marks 
A/=l/3. The inset shows an enlarged view of the kink region 
around h ~ 0.75, illustrating the size dependence. 



D Ei(^ 2 ) 2 [see also Eq. (2)]. We therefore study the ef- 
fect of both frustration J2 and a finite D on the magne- 
tization process of frustrated 5=1 chains. The problem 
of a spin-1 chain with frustration and onsite anisotropy 
D ^ has previously attracted attention as well, how- 
ever, for a slightly different model Hamiltonians than 
our Eq. (2), including biquadratic terms of the form 
(S, • S,) 2 . 15 ' 37 

Our DMRG and ED results for the magnetization pro- 
cess of frustrated 5=1 chains with a finite D are pre- 
sented in Fig. 7 for fixed J2 = 0.6 and A = 1. A 
negative D does not affect the existence of the M=l/3 
plateau state. We have verified the existence of the 1/3 
plateau up to D = — 1. While the width of the plateau 
is Ah « 0.80 at D = and J 2 = 0.6, it increases steadily 
with D < 0. At D = -0.5, we obtain Ah ps 0.95. Thus, 
a negative D stabilizes the M=l/3 plateau state. As an 
example, the magnetization curve M{h) for D = —0.5 
is shown in Fig. 7(a). A positive onsite anisotropy D, 
however, destroys the plateau state which disappears for 
D > 0.18. At larger D, an anomaly around M=2/3 
emerges, as is shown in Fig. 7(d) for the case of D = 0.5. 
Such a behavior at finite D is expected since negative D 
favors the states with a maximal projection 5f on each 
site (such as the uud state on the M=l/3 plateau). On 
the contrary, a positive D plays the role of an easy-plane 
XY anisotropy and thus suppresses the plateau forma- 
tion. 

Let us summarize the main results of this Sec. Ill on 
the magnetization process of frustrated chains with an 
isotropic exchange. First, we have presented numeri- 
cal results for 5=1, 3/2, and 2 showing that a M=l/3 
plateau is realized in the magnetization curve of frus- 
trated spin chains with isotropic exchange, while a clas- 
sical frustrated model does not support this state. Sec- 
ond, we have presented a magnetic phase diagram for 



frustrated 5=1 and 3/2 chains with a special focus on 
the boundaries of the M=l/3 plateau region. A ten- 
dency is found that as 5 > 1/2 grows, the width of the 
plateau becomes more narrow. 

Moreover, kink singularities exist in the magnetization 
curves of spin-1 chains for J2/J1 > 0.25, which is - at 
least close to the saturation field - due to the emer- 
gence of two minima in the dispersion relation of the 
one-magnon excitations above the fully polarized state, 
similar to the case of 5=1/2. 55 This effect further indi- 
cates the presence of an incommensurate region for fields 
h > /ikink above the high-field kink singularity. The in- 
clusion of an onsite anisotropy term in the Hamiltonian 
of frustrated spin-1 chains stabilizes the M=l/3-plateau 
when the onsite anisotropy is negative (D < 0), but the 
plateau is quickly destroyed by a positive D > 0.18. 

The magnetization curves of frustrated spin-3/2 chains 
exhibit a rich behavior with several kink singularities and 
the M=l/3 plateau. Around the plateau and at small 
J2, the magnetization increases in steps of A5 = 25. It 
is beyond the scope of this work to fully map out the 
phase diagram of spin-3/2 chains but interesting results 
are expected to emerge from future studies. 



IV. THE MAGNETIZATION PROCESS OF 
ANISOTROPIC FRUSTRATED SPIN CHAINS IN 
THE EASY- AXIS REGIME A > 1 

In this section we compare our results for frustrated 
spin-5 chains with 5 > 1/2 to a classical frustrated spin 
chain. Beyond a critical exchange anisotropy A C (J2), the 
classical frustrated chain supports a M=l/3 plateau. Us- 
ing linear spin-wave theory summarized in the Appendix, 
we derive several expressions for the phase boundaries of 
the classical model and in particular compare the phase 
boundaries of the plateau region to the cases of 5=1/2,1 
and 3/2 chains. The quantum cases are studied by means 
of ED and DMRG. We find that, as 5 increases, the clas- 
sical phase boundaries are rapidly approached. Devia- 
tions between the quantum cases and the classical case 
are most significant close to the critical frustration for 
plateau formation. The case of 5=1/2 is clearly singled 
out: here, the critical frustration and the phase bound- 
aries are very different from spin S 1,1/2 chains and the 
classical model. Finally, we briefly discuss additional 
phases of frustrated spin-1 chains in the easy-axis regime 
A > 1. 



A. Classical frustrated chain 

We perform MC simulations on chains of up to N = 36 
sites with periodic boundary conditions which allow us 
to determine the ground-state at finite magnetic fields as 
outlined in Sec. II. Knowing the states gives us a handle 
on analytical expressions for their energies. 

The closely related problem of a classical Ising-likc 
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FIG. 9: (Color online) Magnetic phase diagram (field h 
vs frustration J 2) of a frustrated classical chain at A=2. 
The phase boundaries of the classical model are given by 
Eqs. (10a)-(14). They are denoted by straight lines. The 
phases are labeled by capital letters (see Sec. IV B): N is a 
Neel phase, C is a canted phase with periodicity three, DN is 
the double-Neel phase, P the M=l/3 plateau state, CN is a 
canted Neel-like phase, and IC is an incommensurate region. 
Numerical estimates of the phase boundaries of the M=l/3 
plateau region for spin 1/2,1, and 3/2 are denoted by circles, 
squares, stars, respectively. Dashed lines are guides to the eye. 
The dot-dashed line is the saturation field h sat from Eq. (7). 
The dotted line schematically indicates the phase boundary 
between phases CN and IC, which we have not determined 
analytically. 



Heisenberg model on a triangular lattice was studied 
by Miyashita 40 and later on compared to 5=1/2 Ising- 
like Heisenberg model on the same lattice. 64 Similar to 
the Ising limit and spin 1/2, 38 the M=l/3 plateau was 
found and results for the magnetic phase diagram were 
reported. Here we focus on the one-dimensional case 
only, but we allow J\ ^ Ji in our study, while Ref. 40 
concentrated on equal exchange constants on all sides of 
the triangles. 

As an example for our MC calculations, we show the 
magnetization curve for J2 = 0.8 and A = 2 in Fig. 8. 
At least five phases can be distinguished. First, a zero- 
field plateau exists, corresponding to a gap of 1.11. The 
ground state has an up-up-down-down (uudd) structure 
in the S z component, while at small J2, a Neel state (up- 
down) is realized. Larger J2 emphasizes the double-chain 
character of our model. An uudd pattern results in Neel 
states on each of the single chains in the limit of J\ — > 0. 

Next, separated by a first-order transition the system 
enters a region where the magnetization increases linearly 
with h. This phase has periodicity p = 3 and is a precur- 
sor of the plateau state: all spins lie in the same plane, 
one points exactly in the — z direction while the other 
two are canted away from the +z direction by angles 9 
and —9 such that 

C0S fl= - J— ( h +a] . (8) 
A + l \,h + J 2 J w 
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FIG. 10: (Color online) Magnetic phase diagram (field h vs 
exchange anisotropy A) of a frustrated classical chain at J2 = 
0.8 obtained by MC simulations. The phase boundaries of 
the classical model are given by Eqs. (10a)-(14). They are 
denoted by straight lines. Results for the phase boundaries 
of the M=l/3 plateau region for spin 1 and spin 3/2 from 
DMRG calculations are included (squares: spin 1, stars: spin 
3/2). The dashed line is the saturation field h aa ,t from Eq. (7). 
Capital letters (C, DN, and P) refer to phases shown in Fig. 9 
and are discussed in the text. H denotes a helical phase. 65 
The dot-dashed line indicates the transition between phases 
H and C. Dotted lines are phase boundaries that are only 
indicated schematically. 

As h is increased, 9 vanishes continuously, so that the 
transition into the plateau state is second order. The 
M=l/3 plateau exists for Ji + J 2 = 1.8 < h < 6.4108. 
The system then enters a high-field region with a com- 
plicated dependence of M on h. Within the error of our 
simulations the phase transition to this state from the 
plateau state is always second-order. The high-field state 
is probably incommensurate with the lattice period. This 
is illustrated by the chain size dependence of the magne- 
tization curve and its small anomalies shown in the inset 
of Fig. 8. Chains of finite length with PBC cannot ac- 
commodate a truly incommensurate state. This picture 
is consistent with our analytical spin- wave calculations. 
From Eq. (6) we obtain that, for the present value of 
J-2 = 0.8 and h = h sat — + , the fully saturated state 
is instable against fluctuations with cosfc = — J1/AJ2 = 
—0.3125. On the other hand, the plateau state is instable 
against the fluctuations with cos(3fc) = 1 as shown in the 
Appendix. This dependence of k on h emphasizes the in- 
commensurability of the high-field phase. In this work, 
we have not attempted to study the k(h) dependence or 
the nature of this state in further detail. 



B. The magnetic phase diagram of the classical 
model 

We derive analytical expressions for the boundaries of 
the plateau of the classical model as well as of several 
surrounding phases as a function of A and J 2- The results 
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are summarized in the magnetic phase diagrams in Figs. 9 
and 10. 

At low fields, we find three phases: a Neel state at small 
J 2 (region N in Fig. 9), the canted commensurate phase 
C, and the double-Neel state DN for J 2 > 0.64 with an 
uudd pattern. According to Eq. (8), the magnetization 
of the C phase at h — is 

M(h = 0) = 1(2 cosfl - 1) = 1 , (9) 



which vanishes only in the isotropic limit A — > 1. There- 
fore, as a non-trivial consequence of the interplay between 
frustration and exchange anisotropy, we observe a finite 
zero-field magnetization in a model with purely antifer- 
romagnetic interactions. 



The boundaries between the commensurate states 
mentioned above are given by: 



ftN-c = \ {(Ji + Ja)(l - A) + V3v/(Ji + J 2 )(l + A) [A(3Ji - 5J 2 ) - {.h + J 2 )]} , 
dn-c = \ {(Ji + Ja)(l — A) + V3V(Ji + J 2 )(l + A) [A(3J 2 - ,h) - {.h + J 2 )]} . 



r 



(10a) 
(10b) 



For the lower boundaries of the M=l/3 plateau we 
find: 



fc N _ P = 2A(Ji-2J 2 ); 
^c-p = Ji + J 2 ; 
/idn-p = A(2J 2 - Ji) . 



(11a) 
(lib) 
(11c) 



Note that the critical field hc-p does not depend on A. 



Above the plateau, we have identified two phases. At 
small J 2 , a canted Neel state is realized (labeled CN in 
Fig. 9), separated by a first-order transition line 

^n-cn = 2v/(Ji - J 2 )(A - l)[A(Ji + J 2 ) + (Ji - Jzj] 

(12) 

from the gapped Neel state N. The P - CN transition 
is also first-order. It occurs at 



ip-cn 



J 2 + A(Ji + J 2 ) + 2 X /2J 1 J 2 (2 + A 2 ) + J 2 (A 2 - A - 2) + J 2 (A 2 + A - 2)] 



(13) 



r 



For larger J 2 and magnetizations M > 1/3, the situa- 
tion becomes more complicated since incommensurability 
arises, as discussed above. For the range of parameters 
that we study, the P-IC phase transition is always second 
order. Hence, the phase boundary can be derived from 
an analysis of spin-wave instabilities around the plateau 
state (see Appendix), 



Jl + J 2 



'P-IC 



[2A - 1 + V(4A 2 + 4A - 7)] . (14) 



Equations (11a) and (13) result in an expression for the 
miminal frustration J 2 cr i t (A) required for the formation 
of the M=l/3 plateau state in the vicinity of A = 2: 



2,crit 

IT 



i 



2A 2 - 2A X /A(A - 1) 



1 - 2A + 5A 2 



(15) 



This implies J 2 . cr it ~ 0.079 for the classical frustrated 
chain at A = 2. 

To conclude this section, we present a cut through the 
phase diagram at fixed J 2 = 0.8 in Fig. 10, now plotting 
the magnetic field h/S vs the exchange anisotropy A. 



The double-Neel phase DN vanishes at A sa 1.485, and 
below this value of A, a helical phase H emerges at zero 
magnetic field (see, e.g., Ref. 65 for a discussion of heli- 
cal order in the isotropic case). The il/=l/3 plateau only 
exists above a critical value of the exchange anisotropy: 
A sa 1.06. We obtain this result by comparing the en- 
ergy of the helical phase H, which - according to our MC 
simulations - extends up to the plateau at small A, to 
the energy of the plateau phase P. Note that a similar 
plot has been presented in Ref. 40 for the triangular lat- 
tice with equal couplings J on each triangle. A phase 
C has not been reported in that work: at J\ = J 2 , we 
consistently also only find the double-Neel phase below 
the plateau, separated by a first-order transition from the 
plateau [see Fig. 9]. 

We emphasize that we have not yet fully explored all 
regions of the magnetic phase diagrams in Figs. 9 and 10 
of the classical frustrated chain. For instance, an analyti- 
cal expression for the transition between phases CN and 
IC still needs to be found. Moreover, additional phases 
may arise in the limit J 2 3> J\. In this parameter region, 
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FIG. 11: (Color online) Magnetization curves of frustrated 
spin-1 chains with an anisotropic exchange (A = 2) for (a) 
■h = 0, (b) J 2 = 0.2, (c) J 2 = 0.4, and (d) J 2 = 0.8. DMRG 
results (straight lines) are for N — 60 sites, the dashed lines 
are ED results (PBC). The capital letters stand for: Neel 
phase N, canted phase C, double-Neel phase DN. 



FIG. 12: (Color online) Magnetization curves for frustrated 
spin-3/2 chains with an anisotropic exchange (A = 2) for (a) 
J 2 = 0, (b) J 2 = 0.2, (c) J 2 = 0.4, and (d) J 2 = 0.8. DMRG 
results (straight lines) are for N = 60 sites, the dashed lines 
are ED results (PBC). The capital letters stand for: Neel 
phase N, canted phase C, double-Neel phase DN. 



we expect the plateau to vanish as the chains decouple. 



C. M—l/3 plateau in the magnetization process of 
spin-1 and spin-3/2 chains: Comparison to the 
classical case 

We now turn to the magnetization process of frustrated 
quantum spin chains with 5=1 and 3/2 in the easy-axis 
regime A > 1. Numerical results for M(h) are presented 
in Figs. 11 and 12 for J 2 = 0,0.2,0.4,0.8 and A = 2. We 
briefly discuss the magnetization process and emergent 
phases starting with the low- field part. The focus here is 
on to what extent the 5=1 and 3/2 chains resemble the 
phases found for the classical chain. Results on the zero- 
field phase diagram of frustrated spin-1 chains and A > 1 
can be found in Refs. 44 and 46. At a finite A > 1 and 
for spin-1, essentially five phases have been identified: a 
Neel state at small J 2 and A > 1.18, the Haldane phase, 
the Double Haldane phase and the Double-Neel phase for 
A > 1.95 on the large J 2 side (see Ref. 44 and further 
references therein). Our results reported below are in 
agreement with this picture. 

We now concentrate on the example of A=2. First, at 
small J 2-, both the 5=1 and 3/2 chains are in the Neel 
phase N, see for instance Figs. 11(a) and 12(a). This 
phase is gapped and therefore, a zero-field plateau ex- 
ists in the magnetization curves. The Neel phase termi- 
nates at J2,n-C ~ 0.29 for 5=1 and Ja.N-C ~ 0.39 for 
5=3/2. For larger J2, an intermediate region C follows, 
which is gapped in the case of 5=1 [see Fig. 11(c)]. This 
region must contain both the Haldane and the Double- 
Haldane phase, with a transition at J2 ~ 0.8. 44 The com- 
mensurabilty and the existence of a gap in the case of 
5=3/2 in this region still needs to be clarified. The width 



in J2 of this region at zero field becomes narrower as 5 
grows, approaching the boundaries of the phase C of the 
classical model (sec Fig. 9). At J2,c-dn ~ 1-02 for 5=1 
and J2,c-dn ~ 0.94 for 5=3/2, the systems undergoes 
another first-order transition and enters the double-Neel 
phase DN, which is gapped for both values of 5 (see 
Figs. 11(d) and 12(d), J 2 = 0.8). 

At small J 2 and h > 0, the N phase terminates at a 
first-order transition, indicated by macroscopic jumps in 
the magnetization curves of both 5=1 and 3/2 chains [see 
the case of J2 = in Figs. 11(a) and 12(a)]. The critical 
frustrations for the formation of the M=l/3 plateau state 
are J2, C rit ~ 0.16 for 5=1 and J2. C rit ~ 0.12 for 5=3/2. 
Interestingly, an EO region (A5t otal = 2) exists around 
the plateau at small J2 in the case of 5=1 [see Fig. 11(b) 
and (c) for J 2 = 0.2 and J2 = 0.4] and extends far below 
the plateau. For 5=3/2, we find a corresponding region 
with A5^ otal = 3. An example is shown in Fig. 12 (b) for 
J2 = 0.2. We are thus led to conjecture that frustrated 
spin-5 chains generally realize regions with A5j otal = 25. 
We have verified that this is the case for spin 5=2 as well. 

The M=l/3 plateau becomes very broad with increas- 
ing J2, as is evident from the case of J2 = 0.4 depicted 
in Figs. 11(c) and 12(c). Similar to the case of isotropic 
chains, boundary induced steps appear on the plateau. 
We have not determined the upper critical frustrations 
for the plateau formation, but we did verify that the 
plateau disappears for large J2 in both the cases of 5=1 
and 3/2. For magnetizations M > 1/3, we mention that 
kink singularities are again found for J 2 > 0.25, see the 
examples in Figs. 11(c) and 12(c). 

We proceed with a comparison of the critical frustra- 
tions J 2 of the M=l/3 plateau state as a function of 5, 
listed in Table II. For 5 > 1/2, a clear trend towards the 
classical result is seen. Numerical results for the phase 
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Spill S J2,crit/Jl 

1/2 0.45 

1 0.16 
3/2 0.12 

2 0.10 
oo 0.079 



TABLE II: Easy-axis regime, A = 2: Lower critical J2,crit for 
the formation of the M=l/3 plateau for spin 1/2, 1, 3/2, and 
2 (DMRG results) and the classical limit S — > oo. 



boundaries of the plateau state for both 5=1 (squares) 
and 5=3/2 (stars) are included in Figs. 9 (h vs J 2 at 
A = 2) and 10 (h vs A at J 2 = 0.8). It turns out that 
the phase boundaries in Fig. 9 converge rapidly with in- 
creasing 5 towards the classical result. For instance, on 
the scale of Fig. 9, the classical result for the transition 
hc-p from the phase C to the plateau cannot be dis- 
tinguished from the quantum cases. The data points for 
5=1 on the boundary hc-p quantitatively deviate from 
the classical result, but the data points for 5=3/2 are 
already very close to the classical result for J 2 w 0.6. 
Deviations are strongest in both the limit of small and 
large J 2 . 

With respect to Fig. 10 (h vs A at J2 = 0.8), we note 
that sizable deviations between the plateau boundaries 
of the classical chain and the quantum cases only appear 
at small exchange anisotropics close to A = 1. While 
the spin-1 chain realizes the plateau in the isotropic limit 
A = 1, the plateau only opens in M(h) of spin-3/2 chains 
beyond a critical 1 < A crit)3 /2 < 1.02 < A crit , c i ass . 

We can finally conclude that the magnetic phase dia- 
gram of frustrated 5=1 and 3/2 chains in the easy-axis 
regime strongly resembles the phases found in the clas- 
sical limit. Examples are the Neel phase N, the double- 
Neel phase DN, and most importantly, the M=l/3 
plateau state. In the latter case, the similarity is not 
only of qualitative nature, but the phase boundaries of 
the plateau rapidly approach the classical result with in- 
creasing 5. We emphasize our findings on increases of the 
magnetization in steps corresponding to AS'totai = 25 
which may indicate interesting binding effects between 
excitations. 



V. SUMMARY 

In this work we have studied the magnetization pro- 
cess of frustrated spin chains. Our numerical results us- 
ing DMRG and ED techniques reveal a complex and rich 
behavior of the magnetization process of 5=1, 3/2, and 
2 chains. We find jumps, zero- and finite-field plateau 
states, as well as kink singularities. Our primary re- 
sult is the M=l/3 plateau state with broken transla- 
tional invariance and an up-up-down pattern in the spin 
component parallel to the field. The existence of this 



plateau indicates the presence of a gap in the spectrum 
of excitations at finite magnetic fields. We have numer- 
ically estimated the phase boundaries of the plateau for 
5=1 and 3/2 chains with isotropic exchange interactions. 
An additional onsite anisotropy term D(Sf) 2 , necessary 
for the description of many Haldanc materials, stabilizes 
the plateau, when D is negative but rapidly destroys the 
plateau state when D is positive. 

The M=l/3 plateau state can be followed down from 
the Ising-limit by decreasing the ratio A/Ji for fixed 
5. For that reason and for its Neel- type of order, it 
has been classified as a classical state, 24 ' 26 but quantum 
fluctuations seem to stabilize it in the case of isotropic 
exchange: a classical frustrated chain does not show a 
plateau here. We have studied the classical frustrated 
spin chain in detail, which realizes the M=l/3 plateau 
in the easy-axis regime. Our work provides several an- 
alytical results for the magnetic phase diagram of the 
classical frustrated chain including the phase boundaries 
of the M=l/3 plateau for A > 1. In the easy-axis regime, 
the phase diagrams of frustrated 5=1 and 3/2 chains do 
not only qualitatively resemble the classical result, but 
also the phase boundaries of the M=l/3 plateau rapidly 
approach the classical result with increasing 5. The ex- 
treme quantum case of 5=1/2 is singled out, the critical 
frustrations arc different from 5 > 1/2 and the classical 
cases both in the isotropic (A = 1) and the easy-axis 
regime. Our conjecture is that the plateau exists for all 
frustrated spin-5 chains in the easy-axis regime, as our 
results for 5=1, 3/2, 2 and the classical case suggest. 

Our results on the emergence of the M=l/3 plateau 
may be of relevance to the materials mentioned in the 
introduction that are suggested to realize zigzag spin-1 
chains: CaV 2 5 (Ref. 48), NaR(W0 4 ) 2 (Ref. 49), and 
Tl 2 Ru 2 7 (Ref. 50). Moreover, the M=l/3 plateau can 
be expected to emerge in the magnetization process of 
all materials that realize a spin-5 Hcisenberg model on 
a triangular lattice, as is well known for spin 1/2 (see, 
e.g., Ref. 66). Indeed, plateaux at both M=l/3 and 
M=2/3 have recently been reported for the 5=1 bilinear- 
biquadratic Heisenberg model on the triangular lattice. 28 

As by-products to the plateau study, and in particular 
close to the plateau region, we have identified additional 
phases. We have obtained the phase diagram of the frus- 
trated spin-1 chain with isotropic interactions, extending 
the results of other studies. 47 We highlight another in- 
teresting finding: all frustrated spin-5 chains seem to re- 
alize regions where the magnetization increases in steps 
corresponding to A5totai = 25. Extensions of our work 
will comprise a full characterization of different phases 
in terms of correlation functions and a more detailed 
study of excitations. A timely subject is the emergence 
of chirally ordered states in frustrated quantum mag- 
nets, which is investigated both experimentally 48,67 and 
theoretically. 43-45,47 We hope that our work will stimu- 
late future research activities in these directions. 
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FIG. 13: Lowest branch of spin-wave excitations in the uud 
plateau state close to instability. The solid line is for the range 
of parameters studied in Sec. IV, showing instability towards 
k = fluctuations. The dashed line is for a larger J2 with an 
incommensurate instability. 
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APPENDIX A: SPIN WAVES IN THE M-l/3 
PLATEAU STATE 

In this appendix, we outline the linear spin-wave the- 
ory that allows us to determine the phase boundaries of 



the classical model. In this appendix, D = 0. We start 
with the equation of motion for the spin operators S p on 
the pth site of the chain, dS p /dt = i[H, S p ], where H 
is the Hamiltonian Eq. (2) and the brackets denote the 
commutator. Explicitly, we obtain the following, 



dS x 

^ = JpqSq ~ A JpqSq + hS^ 

q q 

dS v 

~~TT = ~~ JpqSq + A Sp JpqSq ~ hS* , (Al) 



where J p ^ q = 3\ for q = p ± 1, J2 for q = p ± 2, and zero 
otherwise. The fluctuations in S z are negligible since 
[H, S*] contains only products of the form Sp SI 1 which 
are of second-order in small fluctuations. Here, we are 
only interested in the classical limit which greatly sim- 
plifies the calculations. Hence, in the following, S p are 
viewed as three-dimensional classical vectors of length S. 

It is convenient to divide the uud plateau structure 
into "unit cells", each containing three chain sites. To 
this end, we replace p by a combined index (m, a), where 
p = 3m + a, m = 0,...,N/3 — 1 is a unit cell num- 
ber, and a = 0,1,2 enumerates the sites within the 
unit cell. Correspondingly, in the linear approxima- 
tion, S^Q = S^, « S and S^ 2 w -S, and (Al) 
becomes a system of linear differential equations for 
S^ a and S^ a . Introducing the plane- wave solution 
Sm V a = Sa' v exp[i(fcp — e t)} and a vector notation for the 
amplitudes S = {S%,S V , Sf, Sf, Sf , Sf ), Eq. (Al) can be 
written as A ■ S = 0, where (7 = e lk ) 



A = 



( ie h 

—h ie 

J17* + J27 2 

- J17* - J2I 2 

- J17* - J2I* 2 



J17 + J27* 2 
-Ji7"^27* 2 
ie h 
—h ie 

- J17* - J21 2 



V J17 + ^27* 2 













-J17* - J27 2 


-J17 - J27* 2 

ic 

-h + 2A( J 1 + J 2 ) 



Jil* + J27 2 \ 


J17 + ^27* 2 


h - 2A(Ji + J 2 ) 

/ 



(A2) 



The spin-wave spectrum e k is defined as the solution to 
the equation 

f(e k ,k) ee dctA = 0. (A3) 

We do not give the explicit cumbersome expression for 
/(e,fc). It turns out, that, as a function of k, /(e, k) 



only depends on cos(3fc) which is consistent with the bro- 
ken translational symmetry of the plateau state. In gen- 
eral, f(e, k) is a third-order polynomial in e 2 , the roots of 
which cannot be written in a closed analytical form. We 
can, however, investigate the minima of the function e k 
that is implicitly given by the equation /(e, k) = 0. The 
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cxtrema of arc defined by 

<k df/dk 
dk ~ df/de 



= 0. 



(A4) 



The plateau state becomes instable towards spin-wave 
fluctuations when e™ m — > 0. After some algebra, we ob- 
tain that the corresponding wave vectors are defined by 
cos(3fc) = 1 or 

cos(3fc) = -4 [hJ 2 + 2 A J 2 ( J : + J 2 ) - J\ - 3 J|)] . (A5) 
4j 2 

For the range of parameters addressed in Sec. IV, the 



incommensurate solution Eq. (A5) is not realized, and 
the critical fluctuations correspond to the first solution. 
As an example, the lowest spin-wave dispersion branch, 
corresponding to the parameters close to the instability, 
is shown in Fig. 13. However, for higher J 2 , the incom- 
mensurate fluctuations may become critical as is shown 
by the dashed line in Fig. 13. A more thorough study of 
that region is beyond the scope of the present paper. 

Substituting cos(3/c) = 1 and = in Eq. (A3), we 
obtain the boundaries of the plateau state, by checking 
the stability against small fluctuations. The solutions 
of Eq. (A6) [see below] yield the second-order transition 
boundaries (lib) and (14). 



(h-Ji- Jif[h 2 + h(Ji + J 2 )(l - 2A) + 2(1 - A)(Ji + J 2 ) 2 f 



(A6) 
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